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Quantum mechanics predicts a number of at first sight counterintuitive phenomena. It is there¬ 
fore a question whether our intuition is the best way to find new experiments. Here we report 
the development of the computer algorithm Melvin which is able to hnd new experimental im¬ 
plementations for the creation and manipulation of complex quantum states. And indeed, the 
discovered experiments extensively use unfamiliar and asymmetric techniques which are challenging 
to understand intuitively. The results range from the first implementation of a high-dimensional 
Greenberger-Horne-Zeilinger (GHZ) state, to a vast variety of experiments for asymmetrically en¬ 
tangled quantum states - a feature that can only exist when both the number of involved parties 
and dimensions is larger than 2. Additionally, new types of high-dimensional transformations are 
found that perform cyclic operations. Melvin autonomously learns from solutions for simpler sys¬ 
tems, which significantly speeds up the discovery rate of more complex experiments. The ability to 
automate the design of a quantum experiment can be applied to many quantum systems and allows 
the physical realization of quantum states previously thought of only on paper. 

PACS numbers: 03.65.Ud, 07.05.Fb, 42.50.Tx 


Quantum mechanics encompasses a wide range of 
counterintuitive phenomena such as teleportation mm, 
quantum interference [3] , quantum erasure [3] , and entan¬ 
glement [5Hin]- Despite our struggle to reconcile them 
with our picture of reality, these phenomena serve as 
building blocks for many exciting and useful quantum 
technologies such as quantum cryptography mi [II] , com¬ 
putation mum, and metrology miiisi- A significant 
challenge arises, however, when we try to combine such 
phenomena in order to perform a complex quantum task. 
Understanding the outcome of even a simple combination 
of these quantum building blocks can be daunting for the 
human intuition. Therefore it is natural to ask: Given a 
certain desired property of a quantum system, what com¬ 
bination of quantum building blocks will be successful in 
achieving it? 

In order to answer this question, we develop a classical 
computer algorithm called Melvin, to which we teach 
how these quantum phenomena work, and subsequently 
assign it a specific problem. The machine then takes 
on the task of hnding and optimizing arrangements of 
quantum building blocks that result in a solution. This 
allows us to uncover experimental methods to create an 
array of new types of entangled states previously thought 
to exist only in theory. In addition, it also allows us to 
address the question of how to manipulate such high¬ 
dimensional quantum states, which is key for their use in 
quantum information systems. 

While searching for these experiments, Melvin en¬ 
larges its own toolbox by identifying useful groups of ele¬ 


ments, leading to a significant speed-up in subsequent 
discoveries. The experiments found by our algorithm 
show a departure from conventional experiments in quan¬ 
tum mechanics in that they rely on highly unfamiliar, but 
perfectly conceivable experimental techniques. This pro¬ 
vides some insight into the kind of out-of-the-box think¬ 
ing that is required for creating such complex quantum 
states. 

Our method aims to create and manipulate general 
complex quantum states for which arbitrary transforma¬ 
tions are not known. The algorithm creates experiments 
using experimental accessible optical components that 
can readily be implemented in the laboratory UTIEH]. 
In addition, our algorithm considers multiple degrees of 
freedom of single quantum systems, and can be extended 
to include nonlinear components and states more com¬ 
plex than single photons. This would allow us to investi¬ 
gate many other interesting quantum phenomena such as 
NOON states [19], induced coherence [20l|2T], quantum 
teleportation of more complex systems [2] or quantum 
metrology [THUS]. A complementary field is computer- 
assisted or automated quantum circuit synthesis (QCS) 
[22H26j . where optimal implementations for quantum al¬ 
gorithms are designed from universal sets of known quan¬ 
tum gates. While very powerful in its own right, the tech¬ 
nique of QCS is used for linear qubit networks and usu¬ 
ally requires fault-tolerant quantum computers for the 
implementation of its results. 

The algorithm - The main goal is to develop an al¬ 
gorithm which finds experimental implementations for 
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Figure 1: Working principle of the algorithm. First, an exper¬ 
iment is created using elements from a basic toolbox. Then, 
the quantum state is calculated, and subsequently its prop¬ 
erties are analyzed. Those properties are compared with a 
number of criteria. If these criteria are not satisfied, the algo¬ 
rithm starts over again. However, if the criteria are satisfied, 
the experiment is simplified and reported, together with all 
relevant information for the user. Useful solutions can be 
stored and used in future experiments, which significantly in¬ 
creases the discovery rate of more complex experiments. The 
orange boxes (toolbox and criteria) are adapted when a differ¬ 
ent type of quantum property is investigated, while the rest 
of the algorithm stays the same. 


quantum states or quantum transformations with in¬ 
teresting properties, see fig.([^. Specific possible input 
states and a toolbox of experimentally known transfor¬ 
mations utilizable by Melvin are defined initially. Using 
the elements from the toolbox, the algorithm assembles 
new experiments by arranging elements randomly. Then, 
from the initial state the resulting quantum state and 
transformation is calculated and its properties are ana¬ 
lyzed. Well-defined criteria that are provided by the user 
decide whether the calculated quantum state has the de¬ 
sired properties. If the quantum states properties satisfy 
the criteria, the experimental configuration is simplified 
and reported to the user. Melvin can store the config¬ 
uration in order to use it as a basic building element in 
subsequent trials. By extending the initial toolbox, it 
is learning from experience, which leads to a significant 
speed-up in discoveries of more complex solutions. 

All quantum states are calculated using symbolic alge¬ 
bra. Every experimental element is a symbolic modifica¬ 
tion of the input state. As an example, a 50/50 symmet¬ 
ric non-polarizing beam splitter for photons is described 
by 


BS['?/’, a, b] = tj} 


a[^] ^ ^(bW+^aH]) 
b[^]^^(a[£]+z.bH]) 


where stands for a symbolic replacement followed by 


a) 3-dimensional 

GHZ state 



b) Schmidt-Rank Vector: (10,6,5) 
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Figure 2: Experimental implementations of high-dimensional 
multipartite entangled quantum states, a) The experimental 
implementation for a 3-dimensional 3-partite GHZ-State. If 
Detector T (Trigger) observes a photon in the state |r) = 
(|0) -I- |1)), then the rest of the quantum state is in a GHZ- 
state, which looks like \ i / j ) = |0,0,0) -I- |1,1,1) -|- |2,2,2) (up 
to local transformations). The parity sorter, as described in 
| 29 |. can sort even and odd 0AM modes. The 3-dimensional 
GHZ-state has a Schmidt-Rank-Vector of (3,3,3) (all compo¬ 
nents are symmetrically entangled with the rest of the state), 
b) A more complex experiment is required for higher-order 
Schmidt-Rank-Vectors. The (10,6,5)-state is one example of 
asymmetrically entangled quantum states. The experiments 
are just two examples of 51 implementations found for creat¬ 
ing a variety of different entangled states. 


a list of substitution rules. £ stands for the orbital angu¬ 
lar momentum (0AM) quantum number of the photon, 
and a and b denote the input paths of the beam split¬ 
ter. For simplicity, all other degrees of freedom (such as 
polarization or frequency) are considered to be the same 
for all photons. For example, for the two-photon state 
Ip = a[3]-b[—3] (a and b represent the path of one photon, 
-1-3 and -3 stand for the 0AM of the photon) the beam 
splitter in path a and b will lead to photon bunching, 
BS[^,a, b] —7> (a[— 3]^ -I- 5[3]^), which is the well-known 
Hong-Ou-Mandel effect [3]. By realizing the calculations 
with symbolic algebra, adding new elements or even new 
degrees of freedom is very easy |27j . Furthermore, it al¬ 
lows easy human-readable intermediate forms, important 
for the examination of solutions and the novel techniques 
found by the algorithm. 

Next we demonstrate the working principles using two 
concrete examples. The demonstrations work in the 
regime of photonic quantum experiments, but the algo¬ 
rithm can readily be adapted to (a combination of) other 
systems such as cold atoms [35] . 

Example 1: High-dimensional multipartite entangle¬ 
ment - The Greenberger-Horne-Zeilinger (GHZ) state is 
the most prominent example for non-classical correla¬ 
tions between more than two involved parties, and has 
led to new understanding of the fundamental properties 
of quantum physics |S]. It has been shown recently that 
its generalization to higher dimensions not only has curi¬ 
ous properties m, but that it is a limiting case of a much 
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richer class of non-classical correlations [HI 133131] ■ Those 
new structures of multipartite high-dimensional entan¬ 
glement are characterized by the Schmidt-Rank Vector 
and give rise to new phenomena that only exist if both 
the number of particles and the number of dimensions are 
larger than two. An example of a state with Schmidt- 
Rank Vector (4,2,2) is the asymmetrically entangled 
state IV'4.2,2) = i(|0,0,0) + |l,0,l) + |2,l,0) + |3,l,l)). 
There, the first particle is 4-dimensionally entangled 
with the other two parties, whereas particle 2 and 3 
are both only two-dimensionally entangled with the rest. 
This Master-Slave-Slave configuration is one of the yet 
unexplored features that only exist in genuine high¬ 
dimensional multipartite entanglement, and will be in¬ 
teresting to study in more detail in future. In order to 
make future experimental investigations possible, we aim 
to find high-dimensional multipartite entangled states in 
photonic systems. 

Here, the initial state is created by a double spon¬ 
taneous parametric down-conversion process (SPDC). 
SPDC is a widespread source for experimental genera¬ 
tion of photon pairs. Multiple SPDC processes can pro¬ 
duce multipartite entanglement, as it is well-known for 
the case of two-dimensional polarization entanglement 
[32l [ 33 ] . However, instead of polarization, we use the 
orbital angular momentum (0AM) of photons [5iH37] . 
which is a discrete high-dimensional degree of freedom 
based on the spatial structure of the photonic wave func¬ 
tion. 

The experiments are generated using a set of basic el¬ 
ements consisting of beam splitters, mirrors, dove prism, 
holograms and OAM-parity sorters P71 E5] . The holo¬ 
grams and the dove prisms have discrete parameters cor¬ 
responding to the 0AM and phase added to the beam, 
respectively. These elements are randomly placed in one 
out of six different paths (four of the paths are inputs 
of the two photon pairs and two are empty to increase 
variability). One arm is used to trigger the tri-partite 
state in the other three arms, which leads to roughly 
10^® possible configurations. At the end, a post-selection 
procedure consisting of the coincidence detection of four 
photons in the first four arms yields the final state. 

We calculate the Schmidt-Rank-Vector of the final 
state and select non-trivial ones (i.e. where there are 
no separable parties). Furthermore, for higher usefulness 
in experiments, we demand that the final state is maxi¬ 
mally entangled in its orbital-angular-momentum. If the 
criteria hold, the experiment is reported. 

Melvin runs for roughly 150 hours (on an Intel Core 
i7 notebook with 2,4 GHz and 24 GB RAM), and finds 
51 experiments for states that are entangled in genuinely 
different ways. Among them, we find the first experimen¬ 
tally realizable scheme of a high-dimensional GHZ-state 
m, a generalization of the well-studied two-dimensional 
GHZ-state (fig. Hh)- Furthermore, we find many ex¬ 
periments for different asymmetrically entangled states 


(such as the |'!/' 4 , 2 , 2 ) explained above). In addition, sev¬ 
eral experiments only differ by continuously tunable com¬ 
ponents (e.g. different holograms or triggers), making it 
possible to explore continuous transitions between states 
of different classes of entanglement. 

The resulting experiments contain interesting novel ex¬ 
perimental techniques previously unknown to the au¬ 
thors. For example, in 50 out of 51 experiments, one 
of the four paths that comes directly from the crystals 
has not been mixed with any other arm (arm D in fig. 
[^, and arm T in |^). The reason is that for double 
SPDC events it is possible that the two photon pairs 
come from the same crystal. Leaving one path unmixed 
leads to erasure of such double-pair emission events in 
four-fold coincidence detection. Interestingly, this im¬ 
mediately introduces asymmetry in the final experimen¬ 
tal configuration. A different novelty is introduced when 
more than 6-dimensional entanglement is created begin¬ 
ning from two three-dimensional entangled pairs. This 
is only possible when the 0AM spectra in two crystals 
are shifted with a hologram and combined in a nontrivial 
way (a preliminary stage of the technique can be seen in 
fig. [^, where the spectrum in arm C is shifted in order 
to reach a 10-dimensional output). In other experiments, 
the normalization of the state has to be adjusted in order 
to get a maximally entangled output. As neutral-density 
filters were not part of the toolbox, Melvin instead used 
beam-splitters as a 50%-filter (for example, fig. [^). 

Now we briefly explain the 3-dimensional GHZ-state 
experiment (fig. [^, details in my- Two indepen¬ 
dent SPDC events in two crystals (which produce 3- 
dimensional entangled pairs) allow for nine different 
states in the four arms. The parity-sorter effectively re¬ 
moves all combinations with opposite OAM-parity from 
two crystals (such as |0, 0, —1, 4-1)), which reduces the 
state to five terms. Detection photon A in the trig¬ 
ger state |T) = |0) -I- |1)) leads to a multipartite 
entangled state where photons C and D reside in a 
three-dimensional space and photon B lives in a two- 
dimensional space HTj. The dimensionality of photon 
B is then increased from 2 to 3 in an intricate combina¬ 
tion of photons A and C. Photon A is shifted by -2 0AM 
quanta and combined with photon C at a beam splitter. 
These photons are then detected in the same mode in one 
BS output, which effectively erases the ’’which-crystal” 
information and entangles the remaining three photons 
into a 3-dimensional GHZ state. 

Example 2: High-dimensional cyclic operations and 
learning - In the second example, we are interested in 
high-dimensional cyclic rotations, which are special cases 
of high-dimensional unitary transformations. A set of 
states is transformed in such a way that the last element 
of the set transforms to the first element (for example, 

11) —>■ |2) —)• |3) —>■ 11) is a 3-cycle). Such transformations 
are required in novel kinds of high-dimensional quantum 
information protocols [3E1133] as well as in the creation 
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Figure 3: Realization of an 8-cyclic rotation using polariza¬ 
tion and OAM |0,y) ...|2,7F) 

| —1, y)). In the experiment, a 4-cycIic rotation for pure OAM 
values is used. Within the 4-cyclic rotation, the parity sorter 
[29] mentioned in the main text is used twice. 

of high-dimensional Bell-states. Here, our input is a set 
of high-dimensional states encoded in different degrees of 
freedom (path, polarization, and OAM). While the cre¬ 
ation and verification of high-dimensional entanglement 
in OAM is well known [Ml 133, the knowledge of how 
to perform arbitrary transformations in this degree-of- 
freedom is still lacking. Thus, finding such transforma¬ 
tions in OAM is very important, as it would enable prac¬ 
tical experiments with high-dimensional quantum states 
and find application in high-dimensional quantum infor¬ 
mation protocols. 

The experiments are generated using a set of basic 
elements that consists of polarizing and non-polarizing 
beam splitters, dove prisms, mirrors, holograms and half¬ 
wave plates. These elements are placed in one of three 
different paths (one path is used as an input, and two 
empty paths are added to increase variability). This leads 
to roughly 10^^ different possible experimental configu¬ 
rations. 

The criterion is based on the largest cycle of the trans¬ 
formation: A number of input states (with different po¬ 
larization (horizontal and vertical), OAM (.^=-10 to -flO) 
and paths) is calculated. Then we search subsets of 
modes that are transformed in a closed cycle, as de¬ 
scribed above, and select the largest closed cycle. Melvin 
was able to find the first experimentally realizable 0AM- 
only 4-cyclic transformation, OAM and polarization 3-, 
6 -, and 8-cyclic rotations and up to 14-cyclic rotations 
using OAM, polarization and path (fig. I&mi). 

Complex problems can be solved more efficiently by 
reusing solutions to simpler problems UniEI]: Whenever 
Melvin finds a solution for a simpler system, it mem¬ 
orizes the experimental configuration as new part of its 
initial toolbox [23 . The novel elements in the toolbox can 
be used to construct the next experimental configuration. 
To compare the effectiveness of learning, we analyze the 
algorithm with and without the ability to increase its own 
set of basic elements. We ran the algorithm for 250 hours, 
and only 3 and 4 instances of 4-cyclic and 6-cyclic rota¬ 
tions were found, respectively. Not a single instance of 
a 3-cyclic and an 8-cyclic rotation was found within 250 
hours. However, using the ability to learn new elements 



Figure 4: Comparison of performance with and without the 
ability to learn (log-scale). Green shows the average time 
required in the case where the algorithm can learn useful 
transformations (the algorithm was executed 10 times with 
the same initial conditions). Black shows the time it requires 
without the ability to learn. The experiments for 3-cyclic and 
8-cyclic transformations were not found (within 250 hours) 
without learning, while experiments for 4-cyclic and 6-cyclic 
rotation were found three and four times in 250 hours, re¬ 
spectively. The errors stand for one standard deviation, cal¬ 
culated from the times it took to hnd the solution. Thus, 
autonomously extending the set of useful transformations im¬ 
proves Melvin s performance, which is crucial for scaling to 
more complex experiments. 

we ran the algorithm 10 times (starting with the initial 
toolbox, i.e. without keeping the learned elements), and 
discovered that the 3- and 6-cyclic rotations were found 
on average within 90 minutes (they were always found 
within 3 hours), and the 4- and 8- cyclic rotations were 
found on average within 3,5 hours (in each of the 10 tri¬ 
als, they were found within 8 hours). Thus the ability 
to learn new elements improves the search by more than 
one order of magnitude, suggesting a mechanism for ex¬ 
periments with a higher complexity (fig. |^. 

Conclusion and Outlook - We have shown how a com¬ 
puter can find new quantum experiments. The large 
number of discoveries reveals a way to investigate new 
families of complex entangled quantum systems in the 
laboratory. Several of these experiments are being built 
at the moment in our labs [13 ng. In contrast to human 
designers of experiments, Melvin does not follow intu¬ 
itive reasoning about the physical system, and therefore 
leads to the utilization of many unfamiliar and uncon¬ 
ventional techniques that are challenging to understand. 
The algorithm can learn from experience (i.e. previous 
successful solutions), which leads to a signihcant speed¬ 
up in discoveries of more complex experiments. 

Melvin can be applied to many other questions about 
the creation and manipulation of quantum systems, such 
as the search for more general high-dimensional trans¬ 
formations with different degrees-of-freedom and for dif¬ 
ferent physical systems such as ultra-cold atoms [28] . 
or for efficient generation of other types of important 
quantum systems such as NOON-states m- In order 
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to improve the efhciency of finding solutions, powerful 
techniques from artificial intelligence research can be ap¬ 
plied, such as evolutionary algorithms [42] (where the 
experiment and the resulting quantum state play the 
role of genotype and phenotype, respectively), reinforce¬ 
ment learning techniques jU 031 03] (by implement¬ 
ing a reward-function depending on the closeness of the 
quantum states properties to the desired properties) or 
entropy-based [3S] and big-data methods [35] (in order 
to find more unexpected solutions). 
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SI) High-dimensional GHZ-entanglement 

A straight-forward generalization of a 3-partite GHZ-state to 3 dimensions looks like 

|V') = ^(|0,0,0) + |1,1,1)+ 12,2,2)). 


( 2 ) 


Local transformations do not change the properties of entanglement, thus a 3-dimensional 3-partite GHZ-state can 
be written as 


= ;^ (|a, b, c) + |a, 5, c) -h |a, 5, c)) , 


(3) 


where x A- x A- x with x = {a, b, c}. 


S2) Spontaneous parametric down-conversion (SPDC) 

The input state of example 1 is a double-emission from SPDG, in the form of 


\'>PsPDc) = N ( ^ |+l?yi,--^ b) + l+^c,--^c)^ 

\i=-DC / 


(4) 


with DC being the highest order of SPDC considered, and A, B, C and D are the path of the photons, and N is a 
normalization constant. We post-select on fourfold coincidences. Such a state can be written with DC=1 gives 


\^spdc,ab+cd) = A^ (|0.4, Ob) + | —1a, 1b) + I 1a, —1b) + |0c, Ob) + | — Ic, 1b) + |lc, — 1b))^ 


(5) 


In the examples, we consider DC=1 up to DC=3 (see chapter S4) to generate the state, and we make sure that up to 
DC=25 higher-order terms dont modify the post-selected output. 

S3) Basis elements from the toolbox 

Here we list the symbolic transformations of all elements used in the toolbox. Here, a or & stand for the path of 
the photon, £ and P stand for the photon’s 0AM and polarization. 


Reflection 


Reflection [i/j, a] = ip ^ 


a[£, a[-£, H] 

a[£, V]^i- a[-£, U] 


( 6 ) 


Non-polarizing symmetric 50/50 beam splitter 


I a[£, P]-t ^ (b[£] -f Reflection[a[^, P]]) 
BS[^, a, b] - V’ ^ I P]^ ^ W] + Reflection[b[f, P]]) 


( 7 ) 
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Polarizing beam splitter 


PBS['^, &,h] = Ip ^ 


a[£,H] h[e,H] 
a[i, V]^i- a[-£, F] 
h[e,H] a[i,H] 
h[£, V]-^i- b[-£, -P] 


( 8 ) 


Half-Wave Plate 


HWP[V', a] = V' 


a[i, H] a[i, F] 
a[i, y] ^ -a[£, H] 


(9) 


0AM hologram 


OAMHolo[i/', a, n] = ■!/) <^ {a[£, P] —)■ a[£ + n, P] 


( 10 ) 


0AM hologram-superposition 


OAMHoloSP[^, a, n] 


iP ■ 


a{l, P] ^ (a[£, P] + a[£ + n, P]) 

v2 


( 11 ) 


Dove prism 


DP['(/',a, n] = tp ^ {a[£, P] —>■ e*"^Reflection[a[£, P]] 


( 12 ) 


0AM Parity Sorter 


a, b] = BS[Reflection[Reflection[DP[Reflection[BS['0, a, b], a], a, l],b], b], a, b] 


(13) 


S4) Experimental Implementations for High-dimensional 3-partite Entanglement 

High-dimensional multipartite entanglement can be characterized by the Schmidt-Rank Vector, introduced in [3]. 
In example 1, we search for such states, and find many examples, as shown in fig.([^ below. All Schmidt-Rank-vectors 
(SRV), their experiments, the required trigger and the resulting quantum state are listed here: 


Found experimental implementations 

for Schmidt-Rank Vectors 


5 2,2 6,2,2 7,2,2 g,2,2 9,2,2 1 0,2,2 
8,3,^ 9^ 10,3,2 

10,3,3 



D „ 

D 

D 


0EF1 



ususa 

EKE 

[jjiaB 

EES 

10,4,3 11,4,3 

12,4,3 

1 4,4,4 

[5,4,4 


7*4,4 

8,4,4 

9,4,4 

10,4,4 11,4,4 

12,4,4 




7AZ 

[8,5,2 

9,5,2 

1 10,5,2 

5,5,3 



|8,5,3 

9,5,3 

10,5,3 

11,5,3 

12,5,3 

5,5,4 


w 


9,5,4 

10,5,4 

11,5,4 

12,5,4 

5,5,5 

6,5,5 

7,5,5 

8,5,5 

9,5,5 

10,5,5 

11,5,5 

12,5,5 



6.6.2I 

7^ 

i 8,6,2 

9,6,2 

10,6,2 

11,6,2 

12,6,2 


6,6,3 

7^,3 


9,6,3 

10,6.3 

11,6,3 

12,6,3 

6,6,4 

7,6,4 

^4 

9,6,4 

10,6,4 

11,6,4 

12,6,4 

6,6,5 

7,6,5 

[Sis 

9,6,5 

10.6^5 

11,6,5 

12,6,5 

6,6,6 

7,6,6 

8,6,6 

9,6,6 

10,6,6 

11,6,6 

12,6,6 



7,7,2 

8,7,2 

9,7,2 

10,7,2 

11,7,2 

12,7,2 

not possible 

DC= 1 

DC= 2 

DC = 3 

7,7,3 

8,7,3 

W 

10,7,3 

11,7,3 

12,7,3 

7,7,4 

8,7,4 

9,7,4 

lOJ^ 

11,7,4 

12,7,4 

7,7,5 

8,7,5 

9,7,5 

10.7^ 

11.7.5 

12,7,5 

7,7,6 

8,7,6 

9,7,6 

10,7,6 

11,7,6 

12,7,6 

7,7,7 

8,7,7 

9,7,7 

10,7,7 

11,7,7 

12.7,7 



[wl 

9,8,2 

10,8,2 

11,8,2 

12,8,2 


8,8,3 

9,8,3 

10,8,3 

11,8,3 

12,8,3 

8,8,4 

9,8,4 

10,8,4 

11,8,4 

12,8,4 

8,8,5 

9,8,5 

10,8,5 

11,8,5 

12,8,5 

8,8,6 

9,8,6 

10,8,6 

11,8,6 

12,8,6 

8,8,7 

9,8,7 

10,8,7 

11,8,7 

12,8,7 

8,8,8 

9,8,8 

10,8,8 

11,8,8 

12,8,8 



9,9,2 

10.9.2 

11,9,2 

12,9,2 


9,9,3 

10,9,3 

11,9,3 

12,9,3 

9,9,4 

10,9,4 

IW 

12,9,4 

9,9,5 

10,9,5 

11,9,5 

12,9,5 

9,9,6 

10,9,6 

11,9,6 

12,9,6 

9,9,7 

10,9,7 

11,9,7 

12,9,7 

9,9,8 

10,9,8 

11,9,8 

12,9,8 

9,9,9 

10,9,9 

11,9,9 

12,9,9 



10.10,21 

11,10,2 

12,10,2 


Figure 5: Found experimental implementations for high-dimensional three-partite entangled states. The list shows experimental 
setups with states of different Schmidt-Rank Vectors. Black cells are not possible m- Cells with green, violet and blue filling 
indicate states that can be generated with 1st, 2nd or 3rd order in the 0AM spectrum (see chapter S2). For white cells, no 
experimental realization has been found yet. 


SRV 

Experimental Setup 

Trigger in A 

resulting quantum state 

DC=1 




(2,2,2) 

"OAMHolo [tA , c, -1] " , "LI [XXX, a, c] " 

|1> +|2> 

-FF2[-1]*FF3[-1]*FF4[0] - FF2 [0] *FF3 [0] *FF4 [1] 

(3,2,2) 

"BS[?/;,b,c]", ’'LI[XXX,c,f]" 

|1> 

-FF2[0] FF3[-1] FF4[0] - FF2[-1] FF3[-1J FF4[1J + 

+ FF2[ll FF3[ll FF4[ll 

(3,3,2) 

"LI [ip ,b, c] " 

1-1) +|0) 

-FF2LlJ*FF3LlJ*FF4L-l] - FF2 [OJ ♦FFS [OJ ♦FF4 [OJ - 

- FF2[ll*FF3[-l1*FF4ril 

(3,3,3) 

"LI ,b, cj " , "Reflection LXXX ,aj " , 

"OAMHolo[XXX,a,-2]", "BS[XXX,a,c]" 

|1> + |0> 

-I FF2[-1] FF3[-3] FF4[-1] - I FF2[0] FF3[-2] FF4[0] + 

+ I FF2[1] FF3[1] FF4[1] 

(4,2,2) 

" OAMHolo [t/j ,b, -1] ", "BS [XXX, a, dj " , 

"BS[XXX,b,d]", "LI[XXX,d,f]" 

|0> 

FF2[-1] FF3[-1] FF4[-1] + FF2 [2] FF3[0] FF4[-1] - 

- FF2[1] FF3[-1] FF4[1] + FF2 [0] FF3[0] FF4[1] 

(4,3,2) 

" OAMH 0 I 0 SP 2 [ip,c, -3] " , "LI [XXX, a, c] " 

|1> +|2) 

-FF2[0J*FF3[0J*FF4[-1J - FF2 [-IJ *FF3 [IJ *FF4 [-1J - 

- FF2[-ll*FF3r-3l*FF4[0l - FF2 [-11 *FF3 [-11 *FF4 [11 

(4,3,3) 

"BS [ip, c,f]", " OAMHolo [XXX, c, 1] " , 

"OAMHolo[XXX,a,2]", "BS [XXX,c,f]", 
"BS[XXX.a,d]", "LI[XXX,c,d]" 

|1> +|0) + 
1-1) 

FF2[1]*FF3[-1]*FF4[-1] + FF2 [1] *FF3 [0] *FF4 [0] + 

+ FF2[0]*FF3[-2]*FF4[2] - FF2 [-1] *FF3 [1] *FF4 [3] 

(5.3,2) 

" OAMHolo [ip,^, -4] " , " BS [XXX, b, c] " , 
"LI[XXX,d,f]", "BS[XXX,a,c]" 

1-1) 

FF2[3]*FF3[-1]*FF4[-1] + FF2 [4] *FF3 [0] *FF4 [-1] + 

+ FF2[1]*FF3[5]*FF4[-1] - FF2 [5] *FF3 [-1] *FF4 [1] + 

+ FF2[-11*FF3r5l*FF4ri1 

(5,3,3) 

"0AMHoloSP2 [i/;, a, 5] " , "LI [XXX,a, c] " 

|0> + 1-1) 

-FF2[11*FF3[lJ*FF4[-11 - FF2 [-1] *FF3 [-6] *FF4 [0] - 

- FF2[1]*FF3[-4]*FF4[0] - FF2 [0] ♦FF3 [0] ♦FF4 [0] - 

- FF2[11*FF3[-11*FF4ri1 

(5,4,2) 

" OAMHolo ['(/!,d,lj " , "OAMHolo LXXX, c,-3J ", 

" LI [XXX, a, d] " , "BS [XXX, b, c] " , 

"LI [XXX.a,c]" 

|0> + |1) 

FF2[0]*FF3[-lJ*FF4[0j + FF2 [-2J ♦FF3 [0] ♦FF4 [0] + 

+ FF2[-4]*FF3[2]*FF4[0] + FF2 [-3] *FF3 [-1] *FF4 [1] - 
- FF2[1]*FF3[3]*FF4[1] 

(5,4,3) 

" OAMHolo [t/j , a, 3] " , "BS [XXX, a, b] " , 

"LI [XXX,a,d]" 

|0> + |1) 

-(FF2L2J*FF3L1J*FF4L-1J) - FF2 L3J *FF3 [OJ *FF4 [OJ - 

- FF2[2]*FF3[-1]*FF4[1] + FF2 [-1] *FF3 [0] *FF4 [2] + 

+ FF2[11*FF3[01*FF4[41 

(5,5,2) 

" OAMHolo [ip,c, -4] " , " BS [XXX, c, d] " , 

"LI [XXX,a,d]" 

|0> + |1) 

-(FF2[-1]*FF3[1]*FF4[-5J) - FF2 [-1J*FF3[-1]*FF4[-3J + 

+ FF2[-1]*FF3[5]*FF4[-1] + FF2 [0] *FF3 [4] *FF4 [0] + 

+ FF2[-11*FF3r3l*FF4ri1 

(6,3,3) 

"0AMHolo['(/!,a,4] " , "Reflection[XXX,b] ", 

"BS [XXX,a,d]", "LI [XXX,a,e]", 
"BS[XXX,a,b]" 

1-1) 

FF2[-3]*FF3[-1]*FF4[-1] + FF2 [-3] ♦FF3 [0] ♦FF4 [0] + 

+ FF2[-3]*FF3[1]*FF4[1] - FF2 [-1] *FF3 [1] *FF4 [3] + 

+ FF2[0]*FF3[-1]*FF4[4] + FF2 [-1] ♦FF3 [-1] *FF4 [5] 

(6,4,3) 

" OAMHolo [ip, a, -4] " , " BS [XXX, a, d] " , 

"LI [XXX,a,c]" 

|0) + |3) 

-(FF2L1J*FF3L0J*FF4L-5J) - FF2 [OJ *FF3 [OJ *FF4 L-4J - 

- FF2[-l]*FF3[0]*FF4[-3] + FF2 [-1] *FF3 [-1] *FF4 [-1] + 

+ FF2[0l*FF.3[-4l*FF4r0l + FF5 [-11 ♦FF.B [11 ♦FF4 [11 
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SRV 

Experimental Setup 

Trigger in A 

resulting quantum state 

(6,4,4) 

"0AMHoloSP2 [-0 ,d, 5] " , "BS [XXX, c, d] " , 
’'LI[XXX,b,d] " 

|0> + |1> 

FF2[-l]*FF3[0]*FF4[-5] - FF2 [-1] *FF3 [6] *FF4 [-1] + 

+ FF2[4]*FF3[-1]*FF4[0] + FF2 [6] *FF3 [1] *FF4 [0] - 
- FF2[0l*FF3[5l*FF4[0l - FF2 [-11 ♦FF3 [41 ♦FF4 [11 

(6,5,2) 

"0AMHoloSP2[0 ,d, 6] " , "BS[XXX,c,d] ", 
’'BS[XXX,a,c]", '’LI[XXX,b,f]’' 

1-1) 

-FF2[-l]*FF3[-lJ*FF4[-7]) - FF2l;-ll *FF3 [OJ *FF4 [-6J + 

+ FF2[l]*FF3[-l]*FF4[-5] + FF2 [-1] *FF3 [-7] *FF4 [-1] + 

+ FF2[-ll*FF3[-6l*FF4[0l + FF2 [-11 *FF3 [-51 *FF4 [ll 

(6,5,3) 

"OAMHolo [0, c, 6] " , "BS [XXX, c , d] " , 

’'LI[XXX,b,c] " 

|-1> + |0> 

FF2[1J*FF3[-5J*FF4[-13 + FF2 [6] *FF3 [0] *FF4 [OJ + 

+ FF2[1]*FF3[-7]*FF4[1] - FF2 [1] *FF3 [1] ♦FF4 [5] - 
- FF2[0l*FF3[0l*FF4[6l - FF2 [ll *FF3 [-11 *FF4 [71 

(6,6,2) 

"LI [0,a,d]", "QAMHolo[XXX,c,3]", 

"BS[XXX,c,d]" 

|-1> + |0> 

FF2LlJ*FF3L-2J*FF4L-l] + FF2 [OJ *FF3 [-3] *FF4 [o] + 

+ FF2[1]*FF3[-4]*FF4[1] - FF2 [1] *FF3 [1] *FF4 [2] - 
- FF2[0l*FF3[0l*FF4[3l - FF2 [ll *FF3 [-11 *FF4 [41 

(7,5,3) 

"0AMHoloSP2[0 , c, 4] " , "0AMHoloSP2[XXX,a,7] ", 
"BS[XXX,c,d]", "LI[XXX,a,c]" 

|1> + |4> 

FF2[-1]*FF3[-3]*FF4[-1] + FF2[-1]*FF3[-8]*FF4[0] + 

+ FF2[1]*FF3[-6]*FF4[0] + FF2 [0] *FF3 [0] *FF4 [0] + 

+ FF2[-1]*FF3[-5]*FF4[1] - FF2 [-1] ♦FF3 [1] ♦FF4 [3] - 
- FF2[-ll*FF3[-ll*FF4[5l 

(7,5,4) 

"OAMHolo [0,d,6] " , "OAMHolo [XXX,b,3] " , 

"BS [XXX, a,b] " , "BS [XXX, c,d] " , 

"LI[XXX,a,d]" 

|1) + |0) 

-FF2[-4]*FF3[1]*FF4[-7J - FF2[-4]*FF3[-1]*FF4[-5] - 
- FF2[-l]*FF3[6]*FF4[-4] - FF2 [1] ♦FF3 [6] ♦FF4 [-2] + 

+ FF2[-4]*FF3[7]*FF4[-1] + FF2 [-3] ♦FF3 [6] ♦FF4 [0] + 

+ FF2[-4l*FF3[5l*FF4[ll 

(7,6,2) 

"0AMHoloSP2[0 , c, -8] " , "BS [XXX,c,d] ", 
"0AMHoloSP2[XXX,c,-5]", "LI [XXX,a,c]" 

|l> + |-4> 

-FF2[0]*FF3[0]*FF4[-9] - FF2 [1] *FF3 [1] *FF4 [-9] - 
- FF2[l]*FF3[-5]*FF4[-8] - FF2 [1] ♦FF3 [-1] *FF4 [-7] + 

+ FF2[1]*FF3[9]*FF4[-1] + FF2 [1] *FF3 [3] *FF4 [0] + 

+ FF2[ll*FF3[7l*FF4[ll 

(7,6,3) 

"OAMHolo[0, c, 5] " , "OAMHolo[XXX,a,-2] ", 
"BS[XXX,c,d]", "BS[XXX,a,d]" 

|1> 

-FF2[-l]*FF3[-l]*FF4[-6] - FF2 [-1J *FF3 [0] *FF4 [-5] - 

- FF2[-l]*FF3[l]*FF4[-4] - FF2 [1] ♦FF3 [-6] *FF4[-3] - 

- FF2[0]*FF3[-6]*FF4[-2] + FF2 [-1] *FF3 [-5] *FF4 [0] + 

+ FF2[-ll*FF3[-4l*FF4[ll 

(8,6,3) 

"0AMHoloSP2[0 ,a,6] " , "BS [XXX,a,b]", 
"0AMHoloSP2[XXX,a,-3]", "LI [XXX,a,d]" 

|0> + |-5> 

FF2[-lJ*FF3LlJ*FF4L-lJ - FF2 L6J *FF3 [OJ *FF4 [OJ + 

+ FF2[-1]*FF3[-1]*FF4[1] - FF2 [5] ♦FF3 [0] ♦FF4 [2] - 
- FF2[7]*FF3[0]*FF4[4] + FF2 [0] ♦FF3 [0] ♦FF4 [6] + 

+ FF2[-ll*FF3[0l*FF4[8l + FF2 [ll *FF3 [01 *FF4 [lOl 

DC=2 




(5,5,3) 

"Reflection[0 ,a] " , "LI[XXX,b,c]", 

"OAMHolo[XXX,a,2]", "BS [XXX,a,c]" 

|-1> + |2> 

(-I)*FF2[-2]*FF3[0]*FF4[-2] + I*FF2 [-1] *FF3 [-1] ♦FF4 [-1] - 

- I*FF2[0]*FF3[0]*FF4[0] - I*FF2 [1] *FF3 [3] *FF4 [1] - 

- I*FF2[2l*FF3[0l*FF4[2l 

(6,4,2) 

"OAMHolo[0,a,3] " , "BS [XXX,a,d] ", 

"OAMHolo[XXX,a,4]", "BS [XXX,a,c]", 
"LI[XXX,c,f]" 

|1) 

-FF2[-2]*FF3[1J*FF4[-1J + FF2 [OJ *FF3 [IJ *FF4 [IJ - 

- FF2[1]*FF3[3]*FF4[2] - FF2 [0] ♦FF3 [3] ♦FF4 [3J - 

- FF2[-1J*FF3[3J*FF4[4J - 

- FF2[-2l*FF3[3l*FF4[5l 

(6,5,4) 

"OAMHolo[0,b,4] " , "OAMHolo[XXX,d,3] ", 

"LI[XXX,b,d]", "OAMHolo[XXX,a,6]", 
"BS[XXX,a,b]" 

1-4)+ |5) 

FF2[2J*FF3[1J*FF4[-6J - FF2 [6J *FF3 [-1J *FF4 [-4J - 

- FF2[8]*FF3[-lJ*FF4[-2] - FF2 [5] ♦FF3 [2] ♦FF4 [1] - 

- FF2[5J*FF3[0J*FF4[3J - FF2 [5J *FF3 [-2J ♦FF4 [5J 

(7,5,2) 

"OAMHolo [0,d, 5] " , "BS [XXX,a,d] " , 
"BS[XXX,a,c]", "LI[XXX,b,f]" 

|1> 

-FF2[-1J*FF3[2J*FF4[-7J - FF2 [-1J *FF3 [IJ *FF4 [-6J - 

- FF2[-lJ*FF3[0j*FF4[-5] + FF2 [1] ♦FF3 [1] ♦FF4 [-4J - 

- FF2[-1J*FF3[-2]*FF4[-3J - FF2 [Ij ♦FF3 [4J ♦FF4 [-1] - 

- FF2[-ll*FF3[4l*FF4[ll 

(8,5,3) 

"OAMHolo[0,b, -5] " , "BS [XXX,b,d]", 

"OAMHolo[XXX,a,-1]", "BS[XXX,a,d]", 
"LI[XXX,d,e]" 

|1) 

FF2[3J*FF3[1J*FF4[-3J - FF2 [5J ♦FF3 [-1J *FF4 [IJ + 

+ FF2[7]*FF3[1]*FF4[1] + FF2 [2] ♦FF3 [-2J *FF4 [5J + 

+ FF2[1]*FF3[-1J*FF4[5J + FF2 [Oj *FF3 [Oj *FF4 [5J + 

+ FF2[-ll*FF3[ll*FF4[5l + FF2 [-21 *FF3 [21 *FF4 [51 

(8,5,4) 

"OAMHoloSP2[0,d,5] ", "LI[XXX,a,d]" 

|-2> + |l> 

-FF2[-2J*FF3[-2J*FF4[-2J - FF2[-1J*FF3[IJ*FF4[-1J - 

- FF2[0]*FF3[-2J*FF4[0J - FF2 [-1] ♦FF3 [-1J *FF4 [Ij - 

- FF2[2]*FF3[-2J*FF4[2J - FF2 [-1] ♦FF3 [2] ♦FF4 [3] - 

- FF2[-ll*FF3[0l*FF4[5l - FF2 [-11 *FF3 [-21 *FF4 [71 

(8,5,5) 

"OAMHolo [0, a, -6] " , "OAMHolo [XXX, a, 2] " , 
"BS[XXX,a,c]" 

|2> 

FF2[2J*FF3[-6J*FF4[-2J + FF2[IJ *FF3[-5J *FF4[-2J + 

+ FF2[0]*FF3[-4j*FF4[-2] + FF2 [-1J *FF3 [-3] *FF4 [-2J - 

- FF2[-2J*FF3[-1]*FF4[-1J - FF2 [-2] *FF3 [0] *FF4 [0] - 

- FF2[-2l*FF3[ll*FF4[ll - FF2 [-21 *FF3 [21 *FF4 [21 

(8,6,4) 

"OAMHolo[0, c, -5] " , "BS [XXX,c,d]", 

"LI[XXX,a,c] " 

|-2> + |l> 

-FF2[-2J*FF3[-2J*FF4[-7J - FF2 [OJ *FF3 [OJ *FF4 [-7J - 

- FF2[2]*FF3[2]*FF4[-7J - FF2 [-1] ♦FF3 [1] ♦FF4 [-6J - 

- FF2[-1J*FF3[-1]*FF4[-4J + FF2 [-1] *FF3 [7] *FF4 [-2J + 

+ FF2[-ll*FF3[5l*FF4[0l + FF2 [-11 *FF3 [31 *FF4 [21 

(8,6,5) 

"LIL0 , d, ej " , "OAMHolo[XXX,b,7J", 

"BS [XXX,b,f]", "LI [XXX,d,f]", 

"BS[XXX,b,e]", "OAMHolo[XXX,a,2]", 
"BS[XXX,a,d]" 

1-1) 

FF2[8J*FF3[-1J*FF4[-1J - FF2 [9J *FF3 [IJ *FF4 [OJ - 

- FF2[7]*FF3[1]*FF4[2] - FF2 [6] ♦FF3 [1] ♦FF4 [3J - 

- FF2[5]*FF3[1]*FF4[4] + FF2 [-2J *FF3 [-2] ♦FF4 [8] + 

+ FF2[0J*FF3[0J*FF4[8J + FF2 [2J *FF3 [2J *FF4 [8J 

(8,7,2) 

"Reflection[0 ,a] " , "0AMHoloSP2[XXX,c,9] ", 

"LI[XXX,d,e]", "LI [XXX,a,c]", 

"0AMHoloSP2 [XXX,b,-3]", "LI [XXX,b,e]" 

1-8)+ |1) 

(-I)*FF2[1J*FF3[1J*FF4[-1J - I*FF2[-1J*FF3[-2J*FF4[1J - 

- I*FF2[1]*FF3[-1J*FF4[1J - I*FF2 [-3] *FF3 [0] *FF4 [1] - 

- I*FF2[-5J*FF3[2J*FF4[1J + FF2[2J*FF3[7J*FF4[lJ + 

+ FF2[0l*FF3[9l*FF4[ll + FF2 [-21 *FF3 [111 *FF4 [ll 

(8,7,3) 

"0AMHoloSP2[0 ,a,6] " , "BS [XXX,a,b]", 

"LI[XXX,a,c] " 

1-5)+ |2) 

FF2[-1J*FF3[1J*FF4[-1J + FF2 [-1J *FF3 [-1J *FF4 [l] - 

- FF2[4]*FF3[-2J*FF4[2J - FF2 [6J *FF3 [Oj *FF4 [2J - 

- FF2[8]*FF3[2]*FF4[2] + FF2 [-2J *FF3 [4J *FF4 [2J + 

+ FF2[0l*FF3[6l*FF4[2l + FF2 [21 *FF3 [81 *FF4 [21 

(8,8,2) 

"Reflection[0,d] " , "OAMHolo[XXX,a,4]", 
"OAMHolo [XXX,c,6]", "LI [XXX,b,c]", 
"BS[XXX,b,d]" 

|2> + |3) 

(-I)*FF2[2J*FF3[-2J*FF4[-8J - I*FF2 [OJ *FF3 [-2J *FF4 [-6J - 

- I*FF2[-2J*FF3[-2]*FF4[-4J + I*FF2 [8] ♦FF3 [-2J *FF4 [-2] + 

+ I*FF2[-1J*FF3[7J*FF4[-1] + I*FF2 [6] *FF3 [-2J ♦FF4 [Oj - 

- I*FF2[ll*FF3[7l*FF4[ll + I*FF2 [41 *FF3 [-21 *FF4 [21 

(9,5,5) 

"0AMHoloSP2[0 , c, 4] " , "BS[XXX,c,d] ", 
"0AMHoloSP2[XXX,a,5] ", "LI[XXX,a,c]" 

|4> + |7> 

FF2[-2J*FF3[-3J*FF4[-1J + FF2[-1J*FF3[-6J*FF4[OJ + 

+ FF2[1]*FF3[-4J*FF4[0J + FF2 [-2] ♦FF3 [-2J ♦FF4 [Oj + 

+ FF2[0]*FF3[0]*FF4[0] + FF2 [2] ♦FF3 [2] ♦FF4 [Oj + 

+ FF2[-2J*FF3[-5]*FF4[1] - FF2 [-2J ♦FF3 [Ij *FF4 [3J - 
- FF2[-2l*FF3[-ll*FF4[5l 

(9,7,3) 

"OAMHolo[0, a, 5] " , "BS [XXX,a,f ] ", 

"OAMHolo[XXX,b,1] ", "BS [XXX,b,f]", 
"BS[XXX,a,c]", "LI[XXX,c,f]" 

|0) + |3) 

-FF2[-3J*FF3[-1J*FF4[-1J + FF2[IJ*FF3[-7J*FF4[OJ + 

+ FF2[-1J*FF3[-5]*FF4[0] + FF2 [-3J *FF3 [-3] *FF4 [0] - 

- I*FF2[4]*FF3[-2J*FF4[0J - I*FF2 [6J ♦FF3 [Oj ♦FF4 [Oj + 

+ I*FF2[-2J*FF3[4J*FF4[0J + I*FF2 [Oj *FF3 [6J *FF4 [Oj - 

- FF2[-3l*FF3[ll*FF4[ll 

(9,7,4) 

"OAMHolo[0,d,3] " , "Reflection[XXX,d] ", 
"OAMHolo[XXX,b,-6]", "BS[XXX,a,b]", 

"LI[XXX,a,d]" 

1-5)+ |4) 

(-I)*FF2[-1J*FF3[-2J*FF4[-5J - I*FF2[-1J*FF3[OJ*FF4[-3J + 
+ I*FF2[4]*FF3[-lJ*FF4[-2] - I*FF2 [-1]♦FF3 [2J♦FF4 [-1] + 

+ I*FF2[6]*FF3[-1J*FF4[0J + I*FF2[8J*FF3[-1]*FF4[2] - 

- I*FF2[-2J*FF3[-1]*FF4[4] - I*FF2 [0] *FF3 [-1J ♦FF4 [6J - 

- I*FF2[2l*FF3[-ll*FF4[8l 
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SRV 

Experimental Setup 

Trigger in A 

resulting quantum state 

(10,5,5) 

" OAMHolo [t/’ ,b, 6] " , "BS [XXX,b, c] " , 
'’LI[XXX,c,e]", "BS[XXX,a,c] " 

|i> 

FF2[2]*FF3[-7]*FF4[-2] + FF2[1]*FF3[-7]*FF4[-1] - 

- FF2[-7]*FF3[1]*FF4[-1] + FF2 [0] *FF3 [-7] *FF4 [0] + 

+ FF2[-1]*FF3[-7]*FF4[1] + FF2 [-8] *FF3 [-2] *FF4 [1] + 

+ FF2[-6]*FF3[0]*FF4[1] + FF2 [-5] *FF3 [1] *FF4 [1] + 

+ FF2r-4l*FF3[2l*FF4[ll + FF2 [-21 *FF3 [-71 *FF4 [21 

(10,6,5) 

" OAMHolo -5] " , " BS [XXX, c, d] " , 

"BS [XXX,b,e]", "BS [XXX,b,f]", 
"BS[XXX,d,e]", "LI[XXX,b,d]" 

|i) 

FF2L-4J*FF3L-lJ*FF4L-l] + FF2[-6]*FF3[ij*FF4[-1] - 

- FF2[2]*FF3[3]*FF4[-1] - FF2 [-1] *FF3 [4] *FF4 [-1] - 

- FF2[0]*FF3[5]*FF4[-1] - FF2 [-2] *FF3 [7] *FF4 [-1] - 

- FF2[-1]*FF3[6]*FF4[1] + FF2 [-1] *FF3 [-2] *FF4 [3] + 

+ FF2[-11*FF3[01*FF4[51 + FF2 [-11 *FF3 [21 *FF4 [71 

(10,6,6) 

" OAMHolo -5] " , " BS [XXX, c, d] " , 

"LI[XXX,b,d]" 

1-2) -til) 

-FF2L-2J*FF3L-7]*FF4[-2J - FF2 [O] *FF3 [-SJ *FF4 L-2J - 

- FF2[2]*FF3[-3]*FF4[-2] + FF2 [-6] *FF3 [-1] *FF4 [-2] + 

+ FF2[-4]*FF3[l]*FF4[-2] - FF2 [-1] *FF3 [-4] *FF4 [-1] - 

- FF2[-1]*FF3[-6]*FF4[1] + FF2 [-1] *FF3 [2] *FF4 [3] + 

+ FF2[-ll*FF3[0l*FF4[5l + FF2 [-11 *FF3 [-21 *FF4 [71 

(10,7,5) 

" 0AMHoloSP2 [i/f, d, 6] " , "BS [XXX, c , d] " , 
"LI[XXX,b,d]" 

1-2) -til) 

FF2[-lJ*FF3LlJ*FF4L-7J + FF2 [-IJ *FF3 [-IJ *FF4 [-5] + 

+ FF2[4l*FF3[-2]*FF4[-2] + FF2 [6] *FF3 [0] *FF4 [-2] + 

+ FF2[8]*FF3[2]*FF4[-2] - FF2 [-2] *FF3 [4] *FF4 [-2] - 

- FF2[0]*FF3[6]*FF4[-2] - FF2 [2] *FF3 [8] *FF4 [-2] - 

- FF2[-ll*FF3[7l*FF4[-ll - FF2 [-11 *FF3 [51 *FF4 [11 

(10,9,2) 

" OAMHolo [t/j , d, 6] " , "BS [XXX, c, d] " , 
"BS[XXX,a,c]", "LI[XXX,b,f]" 

1-1) 

-FF2[-l]*FF3[-2]*FF4[-8] - FF2 [-1J *FF3 [-1] *FF4 [-7J - 

- FF2[-l]*FF3[0]*FF4[-6] + FF2 [1] *FF3 [-1] *FF4 [-5] - 

- FF2[-l]*FF3[2]*FF4[-4] + FF2 [-1] *FF3 [-8] *FF4 [-2] + 

+ FF2[-1]*FF3[-7]*FF4[-1] + FF2 [-1] *FF3 [-6] *FF4 [0] + 

+ FF2[-ll*FF3[-5l*FF4[ll + FF2 [-11 *FF3 [-41 *FF4 [21 

(10,10,2) 

"OAMHolo [t/^ ,d, -9] " , "LI [XXX, a,f ] " , 

"BS [XXX,c,d]", "OAMHolo[XXX,a,2]", 
"BS[XXX,a,d]" 

1-1) 

-FF2[1]*FF3[-2]*FF4[-11] - FF2 [IJ *FF3 [-1] *FF4 [-10] - 

- FF2[l]*FF3[0]*FF4[-9] - FF2 [1] *FF3 [1] *FF4 [-8] - 

- FF2[l]*FF3[2]*FF4[-7] + FF2 [1] *FF3 [-11] *FF4 [-2] + 

+ FF2[1]*FF3[-10]*FF4[-1] + FF2 [1] *FF3 [-9] *FF4 [0] + 

+ FF2[ll*FF3[-7l*FF4[2l - FF2 [-11 *FF3 [-81 *FF4 [31 

(11,7,5) 

" OAMHolo [t/j ,b, 6] " , "QAMHoloSP2 [XXX, d, -5] " , 
"BS[XXX,a,b]", "LI[XXX,a,d]" 

1-1) + |0) 

-FF2[-2J*FF3[0]*FF4[-8] + FF2[-5]*FF3[2]*FF4 [-7] - 

- FF2[0]*FF3[0]*FF4[-6] + FF2 [-5] *FF3 [0] *FF4 [-5] - 

- FF2[2]*FF3[0]*FF4[-4] + FF2 [-5] *FF3 [-2] *FF4 [-3] + 

+ FF2[-8]*FF3[0]*FF4[-2] + FF2 [-5] *FF3 [1] *FF4 [-1] + 

+ FF2[-6]*FF3[0]*FF4[0] + FF2 [-5] *FF3 [-1] *FF4 [1] + 

+ FF2[-41*FF3[01*FF4[21 

DC=3 




(7,4,4) 

"BS[?/;,b,e]", "LI [XXX,b,d] " , 

"BS[XXX,d,e]" 

|3) 

FF2[3]*FF3[-3]*FF4[-3] - I*FF2 [-2] *FF3 [-2] *FF4 [-3] - 

- I*FF2[0]*FF3[0]*FF4[-3] - I*FF2 [2] *FF3 [2] *FF4 [-3] + 
+ FF2[3]*FF3[-1]*FF4[-1] + FF2 [3] *FF3 [1] *FF4 [1] + 

+ FF2[3l*FF3[3l*FF4[3l 

(10,7,4) 

" OAMHolo [t/j , d, 6] " , "BS [XXX, a, d] " , 
"LI[XXX,c,d]" 

|1) + |8) 

FF2[1]*FF3[-3]*FF4[-9J + FF2[1]*FF3[-1]*FF4[-7] + 

+ FF2[l]*FF3[l]*FF4[-5] + FF2 [1] *FF3 [3] *FF4 [-3] + 

+ FF2[l]*FF3[4]*FF4[-2] + FF2 [1] *FF3 [6] *FF4 [0] - 

- FF2[-2]*FF3[-2]*FF4[2] - FF2 [0] *FF3 [0] *FF4 [2] - 

- FF2[2l*FF3[2l*FF4[2l + FF2 [ll *FF3 [81 *FF4 [21 

(11,9,4) 

"LI[?/i,a,d]", "LI[XXX,d,f]", 

"BS[XXX,e,f]", "OAMHolo[XXX,c,5]", 

"BS [XXX,c,d]", "OAMHolo[XXX,c,-4]", 
"BS[XXX,c,e]" 

1-3)-tIO) 

-I*FF2[3]*FF3[6]*FF4[-3] - I*FF2[3]*FF3[8]*FF4[-1] - 

- I*FF2[3]*FF3[10]*FF4[1] + I*FF2 [3] ♦FF3 [1] ♦FF4 [2] - 

- I*FF2[3]*FF3[12]*FF4[3] + I*FF2 [3] *FF3 [3] *FF4 [4] + 

+ I*FF2[2]*FF3[-2]*FF4[5] + I*FF2 [0] *FF3 [0] *FF4 [5] + 

+ I*FF2[-2]*FF3[2]*FF4[5] + I*FF2 [3] *FF3 [5] ♦FF4 [6] + 

+ I*FF2[3l*FF3[7l*FF4[8l 

(12,7,7) 

" OAMHolo [t/^ , a, -6] " , " BS [XXX, a, d] " 

|3) 

FF2[3]*FF3[-3]*FF4[-9] + FF2[2]*FF3[-3]*FF4[-8] + 

+ FF2[l]*FF3[-3]*FF4[-7] + FF2 [0] *FF3 [-3] *FF4 [-6] + 

+ FF2[-l]*FF3[-3]*FF4[-5] + FF2 [-2] *FF3 [-3] ♦FF4 [-4] - 

- FF2[-3]*FF3[-2]*FF4[-2] - FF2[-3]*FF3 [-1]*FF4[-1] - 

- FF2 [-3] *FF3 [0] *FF4 [0] - FF2 [-3] *FF3 [1] *FF4 [1] - 

- FF2[-3l*FF3[2l*FF4[2l - FF2 [-31 *FF3 [31 *FF4 [31 


List 1: Experimental configuration for all found high-dimensional 3-partite entangled states. First column shows the 
Schmidt-Rank Vector. The second column shows the experimental configuration. There, the first element in the list is 
the first in the experiment (for example, the (2,2,2)-states experiment can be calculated as LI[OAMHolo['!/’,c,-l],a,c]). 
The third line is the trigger in path A for the three-photon state in path B, C and D. The last gives the quantum 
state, where FFl, FF2 and FF3 stand for the photon in path B, C and D respectively. 


S5) Example: 3-dimensional GHZ-state (SRV=(3,3,3)) 


The setup for the 3-dimensional GHZ-state in fig. consists of an OAM-Parity sorter, a mirror, a +2 Hologram 
and a beam-splitter. (Note that for simplicity, we do not write normalization-constants) 

After down-conversion (neglecting the double-emissions from one crystal as they will be filtered in the four-fold 
coincidence detection in the end), we can write 
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\^) = (| 0 , 0 ) + | 1 ,- 1 ) + 1 - 1 , 1 ))" = 

= | 0 , 0 , 0 , 0 ) + | 0 , 0 , 1 ,- 1 ) + | 0 , 0 ,- 1 , 1 ) + 

+ | 1 ,- 1 , 0 , 0 ) + | 1 ,- 1 , 1 ,- 1 ) + | 1 ,- 1 ,- 1 , 1 ) + 

+ 1-1,1,0,0)+ 1-1,1,1,-1)+ 1-1,1,-1,1) (14) 

After the OAM-Parity sorter (again, neglecting double-photons in one detector-arm, because they will be filtered by 
four-fold coincidence detection). From here on, only events from the two crystals with the same parity survive: 


\^) = |0,0,0,0) + |1, -1,1, -1) + |1, -1, -1,1) + 1-1,1,1, -1) + 1-1,1, -1,1) (15) 

For a Trigger in A with (|0) + |1)), the state has a SRV of (3,3,2), as one can see in the list above. Photon C and 
D reside in 3 dimensions while photon B lives in a 2-dimensional space. That is because photon B is perfectly anti¬ 
correlated with photon A, which is the 2-dimensional trigger. We want to increase the dimension of photon B to 3. The 
general idea is remove the perfect anti-correlation by mixing the trigger with photon C. At this stage, a BS between 
A and C would lead to Hong-Ou-Mandel interference between the 1st, 3rd and 4th term, which effectively removes 
those terms from the state and does not lead to a 3-dimensional GHZ state. In order to prevent this from happening, 
the photon in A is shifted by -2. This prevents HOM interference between those three terms, and removes the 2nd 
term instead. One additional subtle but significant trick is the usage of the mirror in order to be not vulnerable to 
higher-order terms (without the mirror, the state would become a 2-dimensional GHZ if higher-order modes in SPDG 
are considered). The mirror in arm A leads to: 


= | 0 , 0 , 0 , 0 ) + 1 - 1 ,- 1 , 1 ,- 1 ) + 1 - 1 ,- 1 ,- 1 , 1 ) + | 1 , 1 , 1 ,- 1 ) + | 1 , 1 ,- 1 , 1 ) 


(16) 


And the hologram of -2 in A transforms the state to 


|^) = |-2,0,0,0) + |-3, 


-1,1, -1) + 1-3, -1,-1,1) + 1-1,1,1, -1) + 1-1,1, -1,1) 


(17) 


In the next step, a beam-splitter will be placed between arm A and G. 

IV') = |0,0, -2,0) - |2,0,0,0) + |1, -1, -3, -1) - |3, -1, -1, -1) + |-1, -1, -3,1) - |3, -1,1,1) + 

+ |1,1,-I,-1)-|1,I,-1,-I) + |-1,1,-I,I)-|1,1,1,I) (18) 

The red terms cancel because of destructive interference. This is due to the Hong-Ou-Mandel effect, which occurs if 
the 0AM of two incoming photons from two different arms in a beam splitter are opposite, which leaves us with the 
state 


\xP) = |0,0,-2,0)-|2,0,0,0) + |l,-l,-3,-l) 

+ 1 - 1 , 1 ,- 1 , 1 )-| 1 , 1 , 1 , 1 ) 


|3,-1,-1,-1) + |-1,-1,-3,1)-|3,-1,1,1) + 


(19) 


If we now use the photon A as Trigger for (|0) + |1)), the photons in B, C and D will be in the state: 


IV') =|0,-2,0) + |-l,-3,-l)-|l,l,l) (20) 

This state fulfills the criterion for a high-dimensional GHZ-state as stated above, and with local unitaries can be 
transformed into \'ip) = |0,0,0) + |1,1,1) + |2,2,2). 


6) Cyclic rotations in a high-dimensional space 

Here we list cyclic rotation found by the algorithm, either with 0AM only, or with 0AM and polarization, or with 
0AM and polarization and path. In all examples, V' stands for the initial state and XXX stand for the state after 
the previous element. 
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4-cyclic 0AM rotation 


Experimental configuration: 

BS[V’, a, b] ^ I)F[XXX, b, 1] ^ Reflection[XXX, b] ^ BS[XXX, a, b] 

^ Reflection[XXX, a] ^ BS[XXX, a, b] ^ DP[XXX, b, 1] 

^ Reflection[XXX, b] ^ BS[XXX, a, b] ^ OAMHolo[XXX, a, 1] 


Operation: 


|- 1 )^| 0 )^| 1 )^| 2 )^|- 1 ) 


The number in the ket stands for the 0AM. This experiment has been performed in our laboratories|18j. 


(21) 

( 22 ) 


3-cyclic OAM+Polarisation rotation 


Experimental configuration: 


HWP[i/^, a] ^ Reflection[AAX, a] ^ PBS[AAA, a, c] ^ OAMHolo[AAA, a, 2] 
^ Reflection[AAX, a] ^ PBS [AAA, a, c] ^ BS[AAA, a, b] 

^ DP [AAA, b, 2] ^ Reflection[AAA, b] ^ BS[AAA, a, b] 

^ HWP[AAA, b] ^ BS[AAA, a, b] ^ DP [AAA, b, 2] 

^ Reflection[AAA, b] ^ BS[AAA,a,b] 

Operation: 

|-2,E)^|-2,i/)^|0,E)^|-2,E) 


The number in the ket stands for 


the 0AM, H and V stand for horizontal and vertical polarisation. 


(23) 

(24) 


6-cyclic OAM+Polarisation rotation 


Experimental configuration: 

HWP['(/', a] ^ Reflection [AAA, a] ^ PBS [AAA, a, c] ^ OAMHolo[AAA, a, 2] 
^ Reflection [AAA, a] ^ PBS [AAA, a, c] ^ BS[AAA, a, b] 

^ DP[AAA, b, 2] ^ Reflection [AAA, b] ^ BS[AAA, a, b] 

^ HWP[AAA, b] ^ BS[AAA, a, b] ^ DP[AAA, b, 2] 

^ Reflection [AAA, b] ^ BS[AAA, a, b] 

Operation: 


1-4, H) ^ 1-2, H) ^ [0, E) ^ [2, H) ^ [4, E) ^ [-4, H) 


The number in the ket stands for the 0AM, H and V stand for horizontal and vertical polarisation. 


(25) 

(26) 


8-cyclic OAM+Polarisation rotation 


Experimental configuration: 

PBS[V', a, b] ^ BS[AAA, b, c] ^ DP [AAA, c, 1] ^ Reflection[AAA, c] 

^ BS[AAA, b, c] ^ Reflection [AAA, b] ^ BS[AAA, b, c] 

^ DP [AAA, c, 1] ^ Reflection[AAA, c] ^ BS[AAA, b, c] 

^ OAMHolo[AAA, b, 1] ^ PBS[AAA, a,b] ^ HWP[AAA, a] 


(27) 
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Operation: 


1-1, y) ^ 1-1, if) ^ |o,y) ^ ... ^ \2,H) ^ |-i,y) 


The number in the ket stands for the 0AM, H and V stand for horizontal and vertical polarisation. 


(28) 


14-cyclic OAM+Polarisation+Path rotation 


Experimental configuration: 

Reflection)^:, a] ^ OAMHolo[AAX, a, 2] ^ Reflection[XAA, a] ^ OAMHolo[XAA, a, -2] 
^ PBS[AAA, a,b] ^ HWP[XAA, a] ^ PBS[AXA, a, b] 

^ Reflection[AAX, b] ^ OAMHolo[AAA, a, 2] ^ Reflection [AAA, a] 

^ BS[AAA, a, b] ^ DP[AAA, b, 2] ^ Reflection [AAA, b] 

^ BS[AAA, a, b] 

Operation: 


(29) 


|0,i?,a) ^ |-2,i/,6) ^ |-4,i?,5) ^ |-8,i/,6) ^ |10,y,6) 

^ [-6, H, a) [8, H, a) [6, H, h) [4, H, a) 

—> [0, H, b) —)• [2, V, b) —>■ [2, V, a) —)■ [2, H, a) 

\0,H,a) (30) 

The number in the ket stands for the 0AM, H and V stand for horizontal and vertical polarization, and a and b stand 
for the two different possible paths. 


S7) Learning algorithm 

In the second example involving cyclic operations, the algorithm extends its own set of basic elements autonomously, 
based on the properties of the longest cycle. It saves elements that have large cycles and experiments with non-trivial 
coupling between different degrees-of-freedom. Additionally, elements already learned can be forgotten to improve 
variability and prevent dead-ends, as some of them might even have negative effects on the probability of finding new 
experiments. The decision of which elements are forgotten at which times is purposefully random. Even though it 
would be possible to weight the elements for past usefulness, it would introduce a bias on similar solutions that we 
wanted to prevent. 


S8) Simplification of experiments 

After an experiment is found, it is simplified. For that, three different methods are used. The first one removes 
elements from the experiment and calculates whether the resulting state or transformation is still performed the 
same way. Such simplifications could remove elements in paths that are not accessed. An example, in which it is 
necessary to remove multi elements at the same time, is the following: Four beam splitters after each other form two 
Mach-Zehnder interferometers. Those have no effect if the phases are set correctly, but can only be removed together. 

In the second method, it is tried to replace more complicated elements (such as LI, PBS or DP) by mirrors. This 
works in cases where the only specific modes access the element (for instance, if only vertically polarized photons 
access a PBS). 

A third method tries to simplify the path structure of the experiment by rearranging the paths. For example, if 
two PBSs are used after each other, one output of the second PBS will never be used, thus the second PBS can be 
removed and one path can be removed completely. 

Those three methods are applied iteratively, until no simplification is possible anymore. 



